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A  well-known  textbook  problem  treats  the  mo¬ 
tion  of  a  particle  sliding  frictionlessly  on  the 
surface  of  a  sphere.  An  interesting  variation 
is  to  consider  what  happens  when  kinetic  friction  is 
present. 1  This  problem  can  be  solved  exactly. 

A  free-body  diagram  is  sketched  in  Fig.  1  and  de¬ 
fines  many  of  the  relevant  variables.  The  radial  com¬ 
ponent  of  Newtons  second  law  is 

mgc os  0  —  N  =  mac.  (1) 

Solving  for  N,  which  will  be  needed  to  compute  the 
frictional  force,  we  obtain 

(2) 

r  ) 

where  v(0)  is  the  speed  of  the  object.  Now  we  insert 
this  into  the  tangential  component  of  Newton’s  sec¬ 
ond  law, 

mgs  in  0-  /ulN  =  mat  ,  (3) 

to  obtain 

,  .  „  ^  llv2  dv  dv  d9  //N 

g(smO-ficosO)+J2—  =  —  =  —  — .  (4) 

r  at  aO  at 

Substituting  for  the  angular  speed,  dO! dt  =  vl r,  and 
using  the  identity  IvdvIdO  =  d(v2)ldO  we  get 

^  -  2/jlV2  =  2(sin  0-  /jl  c os  6),  (5) 

dO 

where  the  dimensionless  speed  is  V  =  v  !\frg. 

A  numerical  solution  to  this  differential  equation 
can  be  readily  obtained  using  finite-difference  itera- 
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gc  os 
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Fig.  1.  A  block  of  mass  m  sliding  on  the  surface  of  a 
rough  sphere  of  radius  r,  at  the  instant  it  is  located  at 
angle  6  with  respect  to  the  vertical.  The  three  forces  act¬ 
ing  on  the  object  are  the  normal  force  N,  gravity  mg, 
and  sliding  friction  /jlN,  where  (jl  is  the  coefficient  of 
kinetic  friction.  The  acceleration  has  been  resolved  into 
centripetal  ac  and  tangential  at  components. 


tion  in  a  spreadsheet  for  a  given  initial  speed  and  an¬ 
gle.2  The  equation  also  can  be  solved  analytically  in  a 
calculus-based  course  as  explained  in  Appendix  A. 
This  solution  is  graphed  in  Fig.  2  for  various  values  of 
the  two  parameters  /ul  and  Vq,  where  Vq  is  the  initial 
dimensionless  speed  at  the  top  of  the  sphere.  For  ref¬ 
erence,  curve  1  plots  the  standard  frictionless  example 
of  fi  =  0  and  =  0,  showing  that  the  particle  flies  off 
the  sphere  at  0  =  cos'1  (2/3)  =  48.2°  with  a  dimension¬ 
less  speed  ofV=  A/2/3  =  0.816. 

There  are  two  possible  fates  of  the  object.  It  will  ei¬ 
ther  be  brought  to  rest  on  the  surface  or  it  will  eventu¬ 
ally  lose  contact  with  the  sphere.  The  first  situation 
occurs  if  there  exists  a  value^  of  0  between  0°  and  90° 
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for  which  V=  0,  i.e.,  if  a  plot  of  Kversus  9  hits  the 
horizontal  axis  in  Fig.  2.  Curve  2  shows  a  typical  ex¬ 
ample  of  this  behavior.  On  the  other  hand,  the  object 
flies  off  the  sphere  if  0.  According  to  Eq.  (2),  this 
happens  when  V2  =  cos  0,  i.e.,  if  a  plot  of  Hversus  6 
strays  into  the  shaded  region  of  the  graph  in  Fig.  2. 

(In  particular,  is  constrained  to  be  less  than  1  if  the 
object  is  to  even  begin  on  the  sphere.  This  gives  phys¬ 
ical  significance  to  the  speed  \frg  used  to  normalize 
V.)  For  example,  curve  3  shows  friction  initially  slow¬ 
ing  down  the  particle  but  the  increasing  gradient  of 
the  surface  subsequently  re-accelerating  the  mass 
(which  occurs  in  general  whenever  it  does  not  come  to 
rest  first,  i.e.,  provided  fi  is  not  too  large).  By  starting 
with  just  the  right  speed  for  a  given  value  of  \x ,  the 
particle  can  be  slowed  down  arbitrarily  close  to  zero, 
so  that  the  particle  appears  to  “bounce”  off  the  hori¬ 
zontal  axis  in  Fig.  2.  For  example,  see  curve  4,  dis¬ 
cussed  in  greater  detail  in  Appendix  B. 

The  motion  of  the  block  on  the  surface  of  the  ball 
can  now  be  understood  by  considering  its  trajectory  in 
the  Vwersus-0  parameter  space  of  Fig.  2.  A  curve  be¬ 
gins  at  a  point  on  the  vertical  axis  at  which  V=  V0 
with  0  <  V0  <  1 .  The  plot  proceeds  rightward  until  it 
ends  when  it  contacts  the  horizontal  axis  along  the 
bottom  or  the  gray  region  limiting  its  upward  range, 
whichever  occurs  first.  It  is  an  instructive  exercise  (cf. 
Appendix  B)  to  prove  that  it  is  impossible  for  the  tra¬ 
jectory  to  pass  through  the  intersection  point  of  these 
two  bounding  curves,  i.e.,  the  particle  can  never  reach 
the  equator.  However,  it  can  get  arbitrarily  close  to 
90°  by  skirting  the  shaded  region  all  the  way  along  (cf. 
curve  5  in  Fig.  2). 

In  summary,  a  rich  variety  of  curves  of  v(6)  are 
possible  for  a  point  particle  sliding  on  the  surface  of  a 
rough  sphere.  Generating  and  graphing  these  curves 
can  therefore  prove  a  profitable  method  for  students 
to  learn  how  to  use  a  spreadsheet  in  introductory 
physics.  In  particular,  it  can  help  them  appreciate  that 
while  it  is  helpful  to  try  random  values  of  the  plot  pa¬ 
rameters  (to  generate  curves  2  and  3  in  Fig.  2,  for  in¬ 
stance),  a  more  focused  approach  is  necessary  to  sam¬ 
ple  the  full  range  of  possible  trajectories  (e.g.,  curves  4 
and  5)  of  a  nontrivial  dynamical  system. 

Acknowledgment 

I  thank  John  Mallinckrodt  for  the  definition  of  the 


angle  6  (degrees) 

Fig.  2.  Five  trajectories  of  the  particle  in  V-vs-0  space 
using  the  following  parameters: 


1:  fi  =  0  and  V0  =  0; 

2:  \x  =  0.6  and  V0  =  0.5; 

3:  fJL  =  0.3  and  V0  =  0.6; 

4:  ii  =  1  and  V0  =  0.71944  from  Eq.  (B4); 

5:  n  =  100  and  V0  =  0.9999625  from  Eq.  (B1). 


dimensionless  speed  and  for  suggesting  that  I  investi¬ 
gate  the  limit  as  /ul  — ►  °c  and  — ►  1  simultaneously. 


Appendix  A: 

Analytic  Solution  of  Eq .  (5) 

This  linear,  inhomogeneous,  first-order  ordinary 
differential  equation  can  be  directly  solved  for  V 2  us¬ 
ing  an  integrating  factor.  ^  But  an  alternative  ap¬ 
proach  can  be  used  with  students  who  have  not  yet 
taken  a  differential  equations  course.  One  solution 
can  be  found  by  substituting  the  trial  form 
V2  =  A  cos  9  +  B  sin  9  into  Eq.  (5)  and  separately 
equating  cosine  and  sine  terms  to  find 


Vi 


2 


(4 jul2  -  2)cos  9-6/jl  sin  9 
1  +  4/jl2 


(Al) 


On  the  other  hand,  consider  the  simpler  equation 
obtained  by  setting  the  right-hand  side  of  Eq.  (5)  to 
zero.  This  can  be  rearranged  into  d(V2)l(V2)  = 

2 jmd9  and  both  sides  integrated  to  obtain 


y2  =  Ce2fi0 ? 


(A2) 


where  C  is  an  arbitrary  constant  of  integration.  It  is 
not  hard  to  see  that  V2  +  V2  is  also  a  solution  of  Eq. 
(5) .  (Since  this  sum  contains  one  undetermined  con- 
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stant  and  the  original  differential  equation  is  first 
order,  this  is  in  fact  the  general  solution.  ^)  Finally,  the 
value  of  C  is  obtained  by  fitting  this  general  solution 
to  the  initial  conditions:  The  particle  begins  at  angle 
$0  with  dimensionless  speed  Vq.  In  particular,  sup¬ 
pose  that  0Q  =  0?  he.,  the  object  starts  at  the  top  of  the 
sphere  just  as  in  the  conventional  frictionless  problem. 
Unlike  this  standard  case,  however,  Vq  cannot  be  zero 
because  friction  would  then  cause  the  particle  to 
remain  in  stable  equilibrium  (rather  than  slipping  out 
of  unstable  equilibrium)  at  the  north  pole  of  the 
sphere.^  The  final  solution  therefore  becomes 


U2  = 


(2  -4/jL2)(e2^* * * 8 9  —  cos  6)  —  6/jl  sin  6 

1  +  An2 


+  V$e2»d. 
(A3) 


Appendix  B:  Approaching  the  Equator 


Setting  V=  0  at  6  =  7t/2,  Eq.  (A3)  can  be  solved  for 
V02  to  obtain 


vo2 


4/a2  -  2  + 

1  +  4/a2 


(Bl) 


(In  order  for  this  to  be  positive,  /a  must  be  larger 
than  approximately  0.6034.)  Substitute  this  back 
into  Eq.  (A3)  and  put  6  =  ir/2  —  5  to  obtain 


F2  = 


(An2  -  2)  sin  8  +  6n(e~^  —  cos  8) 

1  +  4/ul2 


(B2) 


By  design,  this  is  zero  at  8  =  0.  However,  it  equals 

(3fi8  -2)8  when  it  is  expanded  to  second  order  in 

8,  i.e.,  V2  is  negative  when  8  is  infinitesimally  small. 
This  means  the  particle  cannot  reach  the  equator:  It 
must  stop  at  some  smaller  angled 

But  how  close  to  the  equator  can  one  get?  In  order 
to  maximize  the  final  angle  9 ,  we  want  d6ld(V2)  0 

at  V=  0.  Consequently  jul  00  according  to  Eq.  (5). 
In  this  limit,  Eq.  (B2)  becomes 


V2  =  sin  8  -  ——  cos  8.  (B3) 

111 

Therefore,  the  object  comes  to  rest  near  7r/2  - 
1.5/ [i.  For  example,  curve  5  shows  what  happens 
when  fi  =  100.  The  particle  actually  stops  at  89.2°, 
in  good  agreement  with  this  prediction. 

This  is  an  unrealistically  large  value  for  jb i.  The  fi¬ 


nal  angle  9  is  maximized  for  a  more  reasonable  value 
of  fi  =  1  when  Vq  is  chosen  so  that  the  trajectory  of 
the  particle  in  Fig.  2  just  grazes  the  horizontal  axis  [so 
that  both  0  and  d{V2)ld0  =  0]  at,  say,  angle  9r 
and  then  speeds  back  up  and  flies  off  the  sphere  at  9. 
The  left-hand  side  of  Eq.  (5)  is  zero,  thus  implying 
that  9r  =  45°,  and  Eq.  (A3)  can  then  be  solved  at  this 
angle  to  deduce  that 

V02  =  0.4^1  +  V2e"nl2j.  (B4) 

Subsequently  the  particle  leaves  the  surface  at  angle  9 
=  69.6°  with  a  dimensionless  speed  of  V  =  0.59, 
plotted  as  curve  4  in  Fig.  2. 
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